Entanglement in XY Spin Chain 



A. R. Its 

Department of Mathematical Sciences, Indiana University-Purdue University Indianapolis, Indianapolis, IN 46202-3216 



o 
o 

(N 

> 

o 

O 

m 
(N 

> 

(N 
O 

o\ 
o 

o 

:^ 



B.-Q. Jin and V. E. Korepin 
C.N. Yang Institute for Theoretical Physics, State University of New York at Stony Brook, Stony Brook, NY 1 1794-3840 

(Dated: February 1, 2008) 

We consider the ground state of the XY model on an infinite chain at zero temperature. Following 
Bennett, Bernstein, Popescu, and Schumacher we use entropy of a sub-system as a measure of 
entanglement. Vidal, Latorre, Rico and Kitaev conjectured that von Neumann entropy of a large 
block of neighboring spins approaches a constant as the size of the block increases. We evaluated 
this limiting entropy as a function of anisotropy and transverse magnetic field. We used the methods 
based on integrable Fredholm operators and Riemann-Hilbert problem. The entropy is singular at 
phase transitions. 



There is an essential interest in quantifying; entan gle- 
ment in various quantum systems |3. |a| 4L|5l |g.[ 7|.la 

lilElElElElElHIilEllillilMlill^ 

tanglement is a fundamental measure of 'quantumness' of 
the system: how much quantum effects we can observe 
and use. It is the primary resource in quantum com- 
putation and quantum information processing 0, Il7| . 
Stable, large scale entanglement is necessary for scalabil- 
ity of quantum computation ji^, 0| • For experimental 
demonstration one can look, for example in 0, E| • The 
XY model in a transverse magnetic was studied from 
the point for view of quantum information in Q , |34l | , 
[53, E^- It was conjectured in JJ that in XY and other 
gapped models the entropy of a block of L neighboring 
spins approaches a constant as L — > 00 (see comment 
|54j). In this paper we evaluated the entropy of a block 
of L neighboring spins in the ground state of XY model 
in the limit L — > 00 analytically, see (|60|l . Hamiltonian 
of XY model can be written as 



n+l 



(l-7)«+i+/i< (1) 



Here 7 is anisotropy parameter (0 < 7 < 1) ; cr^, ct^ cr^ 
are Pauli matrices and /i is a magnetic field. The model 
was solved in ^5], 37J, [s^, 0- Toeplitz determinants 
were used for evaluation of some correlation functions, see 
[4^ and 01- Integrable Fredholm operators were used 
for calculation of other correlations, see [sglliollillli^ l . 

The ground state of the model \GS) is unique. So the 
entropy of the whole infinite ground state is zero, but 
it can be positive for a subsystem [a part of the ground 
state]. We shall calculate the entropy of a block of L 
neighboring spins. We can think that the ground state is 
a binary system \GS) = \AkB). We can call the block 
of L neighboring spins by sub-system A and the rest of 
the ground state by sub-system B. The density matrix 
of the ground state is \GS){GS\. We shall denote it by 
PAB — \GS){GS\. The density matrix of the L neighbor- 
ing spins [subsystem A] is pA ~ TrB{pAB)- Von Neu- 



mann entropy S{pa) of the subsystem A can be repre- 
sented as following: 



S{pa) = -Tr^ (pa In pa), 



(2) 



This entropy defines the dimension of the Hilbert space of 
states of the block of L spins. Majorana operators were 
used in ,3] to describe the entropy (O by the following 
matrix: 



>L = 



Hi Ho 
VHl-i ... 



Ho / 



Here 



g{e) 

-g-\d) 



cos^ — i7sin0 — ft,/2 

and gie) = — . . , . (3) 

' |cos6'-i7sin6l-/i/2| 

One can use an orthogonal matrix V to transform B/, 
to a canonical form: 



1 

-1 



(4) 



The real numbers — 1 < Vm < 1 play an important role. 
We shall call them eigenvalues. The entroiry of a block 
of L neighboring spins was represented in as 

L 

S{pa) - (5) 



with 



H{u) 



.i±i:ini±il-i^lni^. (6) 
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In order to calculate the asymptotic form of the entropy 
let us introduce: 



BL(A) = iA/L-BL, Dl{\)= detBLiX) (7) 



and 



e(x, ly) 



X + V X + ly X — V X — V 
:^ — m — — In — - — . 



(8) 



Here is the identity matrix of dimension 2L. By defi- 
nition, we have H{v) — e(l, v) and 



(9) 



m=l 



In [33| we used Cauchy residue theorem to rewrite for- 
mula jSJ in the following form: 

1 /" d 
S{pa)^ lim —i dAe(l + e,A)— lnI?i(A) . (10) 
e-+o+ 47ri Jyi dA 

Here the contour V in Fig encircles all zeros of Dl{X). 
We also realized that Bl(A) is a block Toeplitz matrix 




FIG. 1: Contours V (smaller one) and V (larger one). Bold 
lines {—oo, —1 — e) and (1 -I- e, oo) are the cuts of integrand 
e(l + e, A). Zeros of Dl{X) (Eq.\^ are located on bold line 
( — 1, 1). The arrow is the direction of the route of integral we 
take and r and R are the radius of circles. ^ 

with the generator ^{z), i.e. 

/Ho n_i ... ni_L \ 

Bi(A)= "° ■ with 

Vhl-i Ho / 



(11) 

, , /At (l-Aiz)(l-A2z-i)\^^^ 
Z-plane with z = x -t- i y 
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FIG. 2: Polygonal line E (direction as labeled) separates the 
complex z plane into the two parts: the part S7+ which lies to 
the left of S, and the part Q- which lies to the right of E. 
Curve H is the unit circle in anti- clockwise direction. Cuts 
Ji,J2 for functions <j){z),w(z) are labeled by bold on line E. 
Definition of the end points of the cuts A... depends on the 
case: Case la: \a = Ai and As = Aj^^, Ac = A2 and A_d = 
Aj^^. Case 16; Xa ~ Ai and Xb = Aj"^, Ac ~ and 
Ad ~ A2. Case 2; Xa ~ Ai and Xb = A2, Ac ~ A^^ and 

Xd^x-\ h 

We fix the branch by requiring that 4>{oo) > 0. We 
use * to denote complex conjugation and S is the unite 
circle shown in Fig. [3 The points Ai, A2 are different 
depending on the Case: 

1. Case la: 2^1-72 <h<2. 



2. Case lb: <h< 2^1-72 

3. Case 2: h>2 . 

In two cases la and 2: 



Ai = 



h- y/h^~ 4(1 - 72) 



A2 — -; -^Ai. 

1 -7 



(13) 



These A are real because /i2 > 4(1 — 72) for both cases. 
In case lb: 



Ai = 



/i -1^/4(1- 72) ~/l2 



A2 = i/xi. 



(14) 



In this case < 4(1 — 72). 

(Note that in the Case 1 the poles of the function ip{z) 
coincide with the points Xa and As, while in the Case 2 
they coincide with the points A^ and Ac-) 
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By virtue of Eq. ifTUI) . our objective becomes the 
asymptotic calculation of the determinant of block 
Toeplitz matrix D]^{X) or, rather, its A -derivative 
^ lnI?i(A). A general asymptotic representation of the 
determinant of a block Toeplitz matrix, which general- 
izes the classical strong Szego theorem to the block ma- 
trix case, was obtained by Widom in ^3 (see also more 
recent work [4^ and references therein). The important 
difference with the scalar case is the non-commutativity 
of the associated Weiner-Hopf factorization. This creates 
serious technical problems. In our work we circumvent 
this obstacle by using an alternative approach to Toeplitz 
determinants suggested by Dcift in [25j. It is based on 
the Riemann-Hilbert technique of the theory of "inte- 
grable integral operators" , which was developed in [2^ , 
[SOj for evaluation of correlation functions of quantum 
completely integrable [exactly solvable] models (see also 
comment ^3)- It turns out that, using the block matrix 
version of j2d| suggested in , one can generalize Deift's 
scheme to the block Toeplitz matrices. In addition, we 
were able to find the explicit Weiner-Hopf factorization 
of the generator ^{z) which eventually made it possible 
to perform an explicit evaluation of the asymptotic of 
the entropy S{pa)- The final result is given in terms of 
elliptic functions and is presented in Eq. 1)65(1 below. In 
what follows we shall outline our calculation providing 
the necessary facts concerning integrable Fredholm op- 
erators. More details, including the evaluation of error 
terms, will be presented in a separate publication. 

Let fj{z) and hj{z), j = 1,2, be 2 x 2 matrix func- 
tions. We introduce the class of integrable operators 
K defined on ^2(2, C^) by the following equations (cf. 

113), 



So we represented Dl{\) as a Fredholm determinant of 
the integral operator K. Define the resolvent operator R 

by 

(I - K){I + R) = I. 

Here / is the identity operator in L2(S,C^). Then we 
have the general equation, 



dX 



IiiDl{\) = -Tv 



which, taking into account that in our case 
d 



dX 



K{z, z') = ~iK{z, z')[h - Hz')r\ 



can be rewritten as 
d 



InDL(A) f tr [i?(z, z)(/2 - *(^))"^] dz. (20) 



dX 



In the formulae above, "Tr" means the trace taking in 
the space L2(S,C^), while "tr" is the 2x2 matrix trace. 
An important general fact is that the resolvent kernel 
satisfies the equation (see e.g. [iof'). 



R{z,z) 



dF^jz ) 
dz 



-H{z). 



(21) 



In this equation, the 4x2 matrix functions F{z) and 
H{z) are determined by the relations. 



F(z) = r+(z)/(z), ze5. 



(22) 



{K X){z)^ <h K{z,z')X{z')dz' for ^ G L2, (15) 



where 



K{z,z') = 



r{z)h{z') 



/(^)- 



/r(-) 



H{z) ^ {Y^y\z)h{z), zeE, (23) 

where the 4x4 matrix function Y^(z) can be found from 
the (unique) solution of the following Riemann-Hilbert 
problem: 

1. Y{z) is analytic for z ^ S. 



Hz) 



hi{z) 
h2{z) 



Let I2 denote the 2x2 identity matrix. Put 

/l(z)=Z^/2, f2{z)^l2 



2TTi 



2iTi 



(16) 



(17) 



hM-z-^K^^ h.iz)^-ll^. (18) 



Then, essentially repeating the arguments of |23 , we have 
the following relation 



DL{X)=Aei{I-K), 



(19) 



2. F(oo) = /4, where 14 denotes the 4x4 identity 
matrix. 

3. Y_{z) = Y+{z).J{z) for z e S where Y+{z) (r_(z)) 
denotes the left (right) boundary value of Y{z) on 
unit circle S (Note: means from inside of the 
unit circle). The 4x4 jump matrix J(z) is defined 
by the equations, 

J(z) =Ii + 2mJ{z)h^{z) 



2/2-$^(z) -z^(/2 - $^(z)) \ 
Z-^{l2-^'^{z)) $^(z) )■ ^^^> 
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Eqs. (|20|l and H21() reduce the original question to the 
asymptotic analysis of the solution Y(z) of the Riemann- 
Hilbert problem (1-3). Our observation is that once again 
we can generalize the arguments of ''2^ to the case of ma- 
trix generator $(z) and arrive to the following asymptotic 
solution of the problem (1-3) ( L 00): 



z^Ul{z)M{z) 
{V^)-\z) 



and 



(y )-^(,)-f iUlr'iz) z^M{z)Vl{z)\ 



(25) 



(26) 



Here 



M(z) = /2-(<i>^)-i(z) 



and U±{z) and V±{z) are 2x2 matrices solving the 
Weiner-Hopf factorization problem : 

(i) <i>{z)^U+{z)U-{z)^V-{z)V+{z), zeS 

(ii) U-{z) and V-{z) {U+{z) and V+{z)) are analytic 
outside (inside) the imit circle S. 

(iii) U-{oo) = V-{oo) = I. 

We can use Eqs. ^ and ^ in Eqs. ^ - ^ and 
obtain the following asymptotic formula: 

±\nDU\)^-j^L + j-l^tT[^{z)]dz, (27) 



*(z) = [U'4z)U^\z) + V^\z)Viiz)] $-i(z), (28) 

as L — > 00 (see also comment j56|'). Here ' means a 
derivative in z variable. 

By explicit calculation, one can find that 



(l-A2)c73$-i(z)a3 = $(z), (T3 



1 

-1 



(29) 



Hence, 



y_(z) = <j3Ul\z)a3 (30) 
V+iz) = a3C/+^(z)a3(l- A2), A^±1, (31) 

and one only needs the explicit expressions for U±{z). 

Our last principal observation is that, for all A out- 
side of a certain discrete subset of the interval [—1,1], 
the solution to the auxiliary Riemann-Hilbert problem 
(i-iii) exists; moreover, the functions U±{z) can be ex- 
pressed in terms of the Jacobi theta-functions. Indeed, 
the auxiliary Riemann-Hilbert problem (i-iii) can be eas- 
ily reduced to a type of the "finite-gap" Riemann-Hilbert 
problems which have already appeared in the analysis 
of the integrable statistical mechanics models (see 



Before we give detail expressions, let us first define some 
basic objects: 

w{z) = ^{z- Ai)(z - A2)(z - A2-i)(z - Ar^), (32) 



1 . A + 1 



2m 



(33) 



where w{z) is analytic on the domain C\ { Ji U J2} shown 
in Fig. El and fixed by the condition: u>(z) z^ as 
z ^ CO. Next we define 



dz 
w(z) 



, c= 2 



dz 
w(z) 



(34) 



Aa 



w{z) ' 



A(z) = / ^ ^ ^ dz, K 



w{z) 



duj{z), (36) 



Points Xa, Ab, Ac, Au and cuts Ji, J2 and curves E and 
S are shown in Fig. |2 We shall also need: 



Ao — lim 

Z — KXj 



A(z)-iln(z-Ai) 



(37) 



Here, the contours of integration for c and S are taken 
along the left side of the cut Ji. The contour of inte- 
gration for T is the segment [As, Ac]. The contours of 
integration for k and in (|37() are taken along the line S 
to the left from A^; also in H37(l . arg(z — Ai) = tt. The 
contours of integration in the integrals A(z) and Lij{z) are 
taking according to the rule: The contour lies entirely in 
the domain il+ (f^_) for z belonging to ri+ {^-)- It also 
worth noticing that ir < 0. 

Now we are ready to introduce the Jacobi theta- 
function, 



3(,S)= J2 



n -\-27Tisn 



(38) 



We remind the following properties of this theta-function 
(see e.g. 

^3(-s) = 03(s), 03(s + l) = 03(s) (39) 
^3(s + t) =e-"'"-2.^%(s) (40) 

6*3 ^n-KmT+ i + = 0, n,meZ (41) 

We also introduce the 2x2 matrix valued function 8(z) 
with the entries, 

en(z) = (z-Ai)-^e^(^) 

es{u;{z) + (3{X)-K+^) 



Uz) 



(TT ' 

2 J 
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ei2(z) = 


-(z-Ai) 2e 






X - 


h {uj{z) - /3(A) 


+ K- 




03 {oj{z) - 


err \ 
2 / 




621(2) = 


-~{z - XiY^er 


A(2) 




X - 


h {u{z) + P{X) 


+ K - 




O3 Hz) - 


¥) 






(z-Ai)-^e'^(^ 


) 




X - 


93 (ojiz) - /3(A) 


- K + 


err \ 

2 y 


03 Hz) + 


2 / 





(42) 



where cr = 1 in Case 1 and cr = in Case 2, and /3(A), 
(x'(z) and K are defined in Eqs. (|32I37|) . The branch of 
(z — Ai)~2 is defined on the 2-plane cut along the part 
of the line E which is to the right of Ai = Aa, and it is 
fixed by the condition arg(z — Ai) = tt, if z — Ai < 0. 

The matrix function & {z) is defined on C\ S. However, 
analyzing the jumps of the integrals uj{z) and A(z) over 
the line S and taking into account the properties H39|) 
and (|40|) of the theta function, one can see that 0(z) 
is actually extended to the analytic function defined on 
C\ { Ji U J2}. Moreover, it satisfies the jump relations 



e+(z) = e-(z)ai 
e+(z) = e_(z)AcriA-i 



z e Ji 
z e J2. 



A = i 
Also note: 



A + 1 \ 








A- 1 / 


CTl = 







611(00) = 

622(00) = 

6i2(^) = 62l(c») = , 



(/3(A) + ^) 

03 (^+^) 

.0^3 (/3(A) -^) 



and 



det6(z) 



(z) dete(oo)W-^. 

Ai 



(43) 
(44) 

(45) 



(46) 

(47) 
(48) 

(49) 



The latter equation follows from the comparison of the 
jumps and singularities of its sides. Finally, we introduce 
the matrix 



Q{z) 



,(z) -</.(z) 
i i 



(50) 



Note that Q{z) diagonalizes original jump matrix $(z): 
$(z) = Q(z)AQ-i(z) (51) 
and Q{z) is analytic on C\ { Ji U J2} and 

Q+iz) = Q^(z)cri, zeJiUJ2- (52) 



We are now ready to present the solution U±{z) of the 
Riemann-Hilbert problem (i-iii). Put 



A^Q{oo)A-^e-^{oo). (53) 



Then, 



U-{z) = AQ{z)AQ-^{z), \z\ > 1 (54) 
U+iz) ^ Q{z)e'\z)A-\ \z\<l. (55) 

Indeed, by virtue of Eq. (|51|l . we only need to be sure 
that U-{z) and U+{z) are analytic for |z| > 1 and |z| < 1 
respectively. From the jump properties of 6(z) and Q{z) 
it follows that U± have no jumps across Ji,2, and hence 
they might have only possible isolated singularities at 
Ai.2,A]/2- The analyticity at these points can be shown 
by observing that the singularities, which the functions 
6(z) and Q{z) do have at the end points of the segments 
Ji^2, are canceled out in the products lf5^ - (|S51) . 

The excluded values of A for which the above construc- 
tion fails are A — ±1 and, in view of Eq. (|49|) . the zeros 
of 03 (/3(A) + ^), i.e. (see (EJ), 



i Am, A^ 



where, 



tanh 771 



l-a 



tttq, m > 0, (56) 



To 



f-^C dz 

.JXb w(z) 
^ — 

r^B dz 
J 'w{z) 



> 0. 



Using explicit formulae H54|l - 155() one can transform our 
basic Eq. I|27() into the form 



^(1- A2) 



tr 



6-(z)-6(z).3 



dz, (57) 



here A 7^ ±1, ±Am- Using the same arguments as for 
Eq. H49|) . one can see that 



1 



cw{z) d(3 



tr 



In 



6-i(z)— 6(z)a3 
dz 



(58) 

This relation allows further simplification of Eq. H27|l . 
Indeed, we have. 



dX 



\nDL{\) 



2A 



tL 



^ 1 



1-A2 

/3(A) + — ) 03 (/3(A) - - 



, (59) 



here A 7^ ±1, ±Am. Taking into account the fact that 
as A ^ 00, Dl{X) — > {-1)^X'^^, we obtain from Eq. 
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the following asymptotic representation for the Toeplitz 
determinant Dl(\): 



D, 



w-|(|)(*^-i)'''>M + t)''»(«^'-t) 



here A lies outside of fixed but arbitrary neighborhoods 
of the points ±1 and ±Am, m > 0. 

It is worth noticing that the asymptotic representa- 
tion for the Toeplitz determinant above shows that, in 
the large L limit, the points A™ (|56|1 are double ze- 
ros of the Dl{^)- This suggests that in the large L 
limit the eigenvalues h'2m and V2m+i from Q merge: 
i^2m, V2m+i — > Am- In turn it indicates the degeneracy of 
the spectrum of the matrix and an appearance of an 
extra symmetry in the large L limit. 

Substituting Eq. (|59|l into the original equation 
Eq. H10|l . and deforming the original contour of integra- 
tion to the contour T as indicated in Fig. ^ we arrive at 
the following expression for the entropy: 



S{pa) 



(60) 



In 



(/3(A) 



2 ) 



(/3(A) 



9i(¥) 



This is a limiting expression as L — > oo. We can 
prove that the corrections in Eq. (|60|l are of order of 

O I^Ap^/VL^ . The asymptotic expression (|6U|I is a the- 
orem, we shall publish a complete proof later. 

The entropy has singularities at phase transitions. 
When r ^ we can use Landen transform (see [32]) to 
get the following estimate of the theta-function for small 
T and pure imaginary s: 



— — s ^irias + O 

IT 



as T 



Now the leading term in the expression for the entropy 
(|^ can be replaced by 



ivr /e^"/^ 



for 



0. 



(61) 



Let us consider two physical situations corresponding to 
small r depending on the case defined on the page 2: 

1. Critical magnetic field: 7 7^ and ft, ^ 2. 

This is included in our Case la and Case 2, when 
h > 2-\/l — 7^. As /i — > 2 the end points of the cuts 
As ^ Ac, so T given by Eq. H34(l simplifies and we 
obtain from Eq. 1)61(1 that the entropy is: 



S{pa) 



iln|2 



- h\ iln47, 
77^0 



(62) 



correction is 0(|2 — /i| In^ |2 — /i|). This limit agrees 
with predictions of conformal approach |3, (3- 
The first term in the right hand side of (|62|l can be 
represented as (1/6) In^, this confirms a conjecture 
ofjl^.^ The correlation length ^ was evaluated in 

2. An approach to XX model: 7 — > and h < 2: It 
is included in Case lb, when < h < 2\J\ — 7^. 
Now As ^ Ac and A^ Ad, we can calculate t 
explicitly. The entropy becomes: 

S\pa) = -\ 1117 + iln(4 - ft') + i ln2, 

7^0 ft < 2 (63) 

correction is 0(7 In' 7). This agrees with js^. 

It is interesting to compare this critical behavior to 
Lipkin-Meshkov-Glick model. It is similar to XY model 
but each pair of spins interact with equal force, one can 
say that it is a model on a complete graph. The crit- 
ical behavior in Lipkin-Meshkov-Glick was described in 
[4^ . it is similar to XY but actual critical exponents are 
different. 

Note. We can integrate over the original contour F' 
of Fig. n after substituting Eq. ((59() into the equation 
Eq. (|10|l . This will give an alternative representation for 
the entropy S{pa) in terms of an infinite series: 



S{pA)^2j2H{X,n)^ J2 (1 + A„01n- 



m=0 



+ A, 



here the numbers A„i are defined in Eq. H56() . This rep- 
resentation is similar to Eq. (jsj. 

Remark. These numbers Am satisfy an estimate: 

|A,„+i - Xm\ < 47rTo with tq = -ir. 

This means that (A,„+i — Xm) — > as r ^ for every 
TO. This is useful for understanding of large L limit of 
the XX case corresponding to 7 — s- 0, as considered in 
[3^ . The estimate explains why in the XX case the 
singularities of the logarithmic derivative of the Toeplitz 
determinant dliiDL(X)/dX form a cut along the interval 
[—1,1], while in the XY case it has a discrete set of poles 
at points ±A„i of Eq. (|56|l . 

It would be interesting to generalize our approach to 
the a new class of quantum spin chains introduced re- 
cently b y J . Keating and F. Mezzadri, while study matrix 
models 



SUMMARY 

Our main result is the theorem that the expression for 
the limit (as L ^ 00) of the entropy of a block of L 
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neighboring spins on the infinite lattice is given by for- 
mula H60|) . We are preparing a larger file with all the 
details of the proof. We can change variables and repre- 
sent it in the form: 



S{pa) = 



(65) 



2 Jo 




We remind that a = for Case 2 and cr = 1 in Case 1, see 
page 2. S{pa) implicitly depends on 7 and h introduced 
in 1^ by means of r defined in Eq. (|34(l , this r is a ratio of 
periods for the theta function 6*3 (z), see (|38|l . The 7 and 
h define A^, Ab, Ac and Ad by means of Eqs. ifT^ . 
Fig- El and its caption. Finally Eq. and Eqs. or 
(El define S{pa) ■ 



APPENDIX 

After our paper appeared in quant-ph, I.Peschel IsJl 
simplified our expression for the entropy in the Cases la 
and 2. He used the approach of He showed that 
in these cases our formula (|t)4|) is equivalent to formula 
(4.33) of "9^. Moreover, I. Peschel was able to sum it up 
into the following expressions for the entropy. 



In 



16/c' 

in Case la, and 



1 

2 



AI{k)I{k') 



In 2, 



S= — 
12 



in Case 2 



In 



16 



+ (A:^ - k'^) 



,2^4/(fc)/(fc') 



(66) 



(fc2fc'2) 

Here, I{k) denotes the complete elliptic 



integral of the first kind, k' = ^1 — k-^, and 

V(V2)2+7^-l / 7, Case la 
7 / ^(V2)2+72-l, Case 2 



k = 



(67) 



In our work, we have shown, in particular, that equa- 
tion H64|l is valied in Case lb as well. Therefore, we can 
apply the summation procedure of ^5l]| and obtain that 
in Case lb 



S = 



6 



In 



/ fc2 



Vl6fc' 



2 



AI{k)I{k') 



In 2 



with k' = Vl — fc^, and 



k = 



'1- (V2)2-72 
1 - (h/2r 



(68) 
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